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A MET HOD AND APPARATUS FOR IMAGING USING S HEAR WAVES 

The present invention relates to imaging methods and 
apparatuses using shear waves . 

More particularly, the invention relates to an 
5 imaging method using shear waves for observing a 

diffusing viscoelastic medium containing particles that 
reflect ultrasound compression waves, said method 
comprising: 

a) an excitation step during which an elastic shear 
10 wave is generated in the viscoelastic medium; 

b) an observation step during which the propagation 
of the shear wave is observed simultaneously at a 
multitude of points in an observation field in the 
viscoelastic medium, this observation step comprising the 

15 following substeps : 

bl) causing an array of transducers that are 
controlled independently of one another to emit into the 
viscoelastic medium a succession of unfocused ultrasound 
compression wave shots at a rate of at least 500 shots 

20 per second; and 

b2) causing sound signals received from the 
viscoelastic medium to be detected and recorded in real 
time, said signals comprising the echoes generated by the 
unfocused ultrasound compression wave interacting with 

25 the reflecting particles in said viscoelastic medium; and 

c) at least one processing step during which: 

cl) the sound signals received successively 
from the viscoelastic medium during substep b2) are 
processed in order to determine successive propagation 
30 images of the shear wave; and 

c2) at least one movement parameter of the 
viscoelastic medium is determined at different points of 
the observation field. 

This produces a "motion picture" clearly 
35 illustrating the propagation of the shear wave through 
the viscoelastic medium, which can make it possible to 
perform qualitative and/or quantitative analysis in order 
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to identify zones having hardness that differs from the 
hardness of the remainder of the viscoelastic medium, or 
zones having relaxation time that differs from the 
relaxation time of the remainder of the viscoelastic 
5 medium . 

Document WO-A- 00/55616 describes an example of such 
a method, in which shear waves are generated at the 
surface of the viscoelastic medium. That method gives 
full satisfaction in particular when imaging zones 
10 -situated relatively close to the surface of the 

viscoelastic medium. However that known method does not 
enable certain zones to be observed in the viscoelastic 
medium, and in particular: 

- zones that are sufficiently deep to be unreachable 
15 by shear waves generated at the surface (shear waves 

attenuate quickly) ; and 

- shadow zones that are masked by obstacles (in 
particular portions of a patient's skeleton or liquid 
zones such as liquid cysts) which impede the propagation 

2 0 of shear waves. 

In addition, if the observation field is partially 
in a shadow zone, it can be necessary to move the shear 
wave generator device during an observation, which is 
tedious for the user. 

25 Finally, the shear wave generator device is 

relatively heavy and complicates the apparatus. 

A particular object of the present invention is to 
mitigate those drawbacks. 

To this end, according to the invention, a method of 

30 the kind in question is characterized in that during 

excitation step a) the elastic shear wave is caused to be 
generated by causing at least one focused ultrasound wave 
to be emitted into the viscoelastic medium by said array 
of transducers, the focusing and the timing of said 

35 focused ultrasound wave, and the timing of said unfocused 
ultrasound wave being adapted so that at least some of 
saiid unfocused ultrasound waves penetrate into the 
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observation field while the shear wave is propagating in 
the observation field- 

Thus, the same array of transducers can be used both 
for generating the elastic shear wave in selected manner 
5 in the observation field, and for subsequently, observing 
said propagation by virtue of the fact that the imaging 
apparatus is adapted to generate either focused 
ultrasound waves enabling the elastic shear waves to be 
generated or unfocused ultrasound waves enabling the 
10 propagation of the shear wave to be observed, and 
suitably selecting: 

- the timing of the various emissions; and 

- the point (s) on which the focused ultrasound wave 
is/are focused. 

15 The imaging method of the invention is thus easy to 

implement for a user using apparatus that is relatively 
simple and lightweight. The invention is thus of very 
low cost compared with competing techniques such as 
magnetic resonance imaging (MRI) , and where appropriate 

2 0 it makes it possible to establish outpatient imaging 
systems that can be used for preoperative imaging, 
postoperative imaging, and even for imaging while an 
operation is in progress. 

By way of example, in medical applications, the 

25 method of the invention can make it possible for 

cancerous zones within the tissues of a patient to be 
identified effectively. Shear waves propagate through 
cancerous zones in very different manner than through 
adjacent zones. This identification can bei performed 

30 much more easily than by conventional observation using 
simple ultrasound echography, since the propagation of 
shear waves is a function of the shear modulus of the 
medium, which is itself highly variable between zones of 
healthy tissue and a zone of cancerous tissue: typically, 

35 the shear modulus varies over a ratio of 1 to 30 between 
a healthy zone and a cancerous zone, whereas the 
compression modulus, which governs the propagation of 
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compression soundwaves as used in ultrasound echography 
varies only by about 5% between healthy tissue and 
cancerous tissue. 

Similarly, it is thus possible to identify zones of 
5 necrosis within tissue, for example tumor zones that have 
been subjected to ultrasound hyperthermia treatment, in 
particular for the purpose of evaluating the 
effectiveness of the hyperthermia treatment. 

Another possible application of the invention 

10 relates to quantitatively evaluating the degree of 

fibrosis of the liver, which is an important parameter in 
liver disease, in particular hepatitis C. 

It shbuld be observed that the invention makes it 
possible to generate the shear wave, and to observe its 

15 propagation, including through a zone of liquid or 

through a bone barrier (skull, rib cage, etc.) completely 
or partially masking the observation field, since it is 
possible to focus ultrasound waves through such barriers 
(see in particular document WO-A-02/32316 or French 

20 patent application No. 02/10682 of August 28, 2002) . 

In preferred implementations of the method of the 
invention, it is optionally possible also to have 
recourse to one or more of the following dispositions: 
- during substep b2) , in order to determine said 

25 movement parameter, a plurality of successive propagation 
images (e.g. by correlation, Doppler, etc.) are compared 
with a common reference image of the viscoelastic medium, 
the reference image being determined by firing at least 
one unfocused ultrasound compression wave into said 

3 0 viscoelastic medium and then detecting and recording 

echoes generated by said unfocused ultrasound compression 
wave on interacting with the reflecting particles in the 
viscoelastic medium (thus improving the accuracy with 
which the movement parameter is measured (e.g. 

35 displacement) in the viscoelastic medium, in particular 
for movement of small amplitude (typically less than 



5 



30 micrometers (/im) with the soundwave excitation 
technique used herein) ) ; 

- step a) is preceded by an initial observation step 
aO) during which at least one unfocused ultrasound 

5 compression wave is fired and then echoes generated by- 
said unfocused ultrasound compression wave interacting 
with the reflecting particles in the viscoelastic medium 
are detected and recorded, said echoes corresponding 
(directly or indirectly) to an initial image of the 
10 viscoelastic medium, and during substep b2) , said initial 
image constitutes said reference image for processing at 
least some of the successive displacement images; 

- during initial observation step aO) , a plurality 
of unfocused ultrasound compression waves are fired in 

15 succession and then echoes generated by each unfocused 
ultrasound compression wave interacting with the 
reflecting particles of the viscoelastic medium are 
detected and recorded, said echoes corresponding 
(directly or indirectly) to a plurality of successive 

20 images of the viscoelastic medium, and said initial image 
of the viscoelastic medium is determined by combining 
said successive images; 

- said movement parameter is a displacement of the 
viscoelastic medium; 

25 - the focused ultrasound wave emitted during 

excitation step a) presents a frequency f_ lying in the 
range 0.5 megahertz (MHz) to 15 MHz, and is emitted for a 
duration of k/f seconds, where k_ is an integer lying in 
the range 50 to 5000 and f_ is expressed in hertz (Hz) ; 

30 - the focused ultrasound wave emitted during 

excitation step a) presents a frequency lying in the 
range 0.5 MHz to 15 MHz and is emitted during a 
succession of emission periods separated by rest periods, 
the emission periods following one another at a rate 

35 lying in the range 10 to 1000 emissions per second; 

- the focused ultrasound wave emitted during 
excitation step a) is a linear combination (in particular 
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a sum) of two monochromatic signals having respective 
frequencies fl and f2 such that 20 Hz < |fl - f2| < 
1000 Hz; 

- the focused ultrasound wave emitted during 
5 excitation step a) is focused simultaneously on a 

plurality of points; 

- image processing step c) is followed (immediately 
or otherwise) by a mapping step d) during which, on the 
basis of variation in the movement parameter over time, 

10 at least one shear wave propagation parameter is 

calculated at at least some points of the observation 
field in order to determine a map of said propagation 
parameter in the observation field; 

- the shear wave propagation parameter which is 

15 calculated during mapping step d) is selected from shear . 
wave speed, shear modulus. Young's modulus, shear wave 
attenuation, shear elasticity, shear viscosity, and 
mechanical relaxation time; and 

- steps a) to d) are repeated successively while 
20 emitting different focused ultrasound waves during 

successive excitation step a) , and then combining the 
maps obtained during the successive mapping step d) in 
order to calculate a combination map of the observation 
field. 

25 Furthermore, the invention also provides an imaging 

apparatus for implementing a method according to any 
preceding claim using shear waves to observe a diffusing 
viscoelastic medium containing particles that reflect 
ultrasound compression waves, the apparatus comprising an 

30 array of transducers that are controlled independently of 
one another by at least one electronic central unit 
adapted: 

- to cause at least one elastic shear wave to be 
generated in the viscoelastic medium; 

35 - to observe the propagation of the shear wave 

simultaneously at a multitude of points in an observation 
field in the viscoelastic medium by causing said array of 
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transducers to emit into the viscoelastic medium a 
succession of unfocused ultrasound compression wave shots 
at a rate of at least 500 shots per second^ then causing 
said array of transducers to detect in real time and 
5 record in real time sound signals received from the 
viscoelastic medium, the sound signals comprising the 
echoes generated by the unfocused ultrasound compression 
wave interacting with the reflecting particles of said 
viscoelastic medium; and 

10 - processing the successive sound signals received 

from the viscoelastic medium to determine successive 
propagation images of the shear wave, and then 
determining at least one movement parameter of the 
viscoelastic medium at different points of the 

15 observation field; the apparatus being characterized in 
that the electronic central unit is adapted to cause the 
elastic shear wave to be generated by causing at least 
one focused ultrasound wave to be emitted into the 
viscoelastic medium by said array of transducers, the 

20 focusing and the timing of said focused ultrasound. wave, 
and the timing of said unfocused ultrasound wave being 
adapted so that said unfocused ultrasound waves reach the 
observation field during the propagation of the shear 
wave through the observation field, 

25 Other characteristics and advantages of the 

invention appear from the following description of an 
embodiment thereof, given by way of non-limiting example 
and with reference to the accompanying drawing. 

In the drawing. Figure 1 is a diagrammatic view of a 

30 shear -wave imaging device in an embodiment of the 
invention . 

The imaging device 1 shown in Figure 1 is for 
studying the propagation of elastic shear waves in a 
viscoelastic medium 2 that diffuses ultrasound waves in 
35 compression, and that may be constituted, for example: 
- by an inert body, in particular for quality 
control in industrial applications; or 
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- a living body, for example a portion of the body 
of a patient, in medical applications. 

By way of example, these movements are tracked by 
means of a microcomputer 4 (comprising at least an input 
5 interface 4a such as a keyboard, etc., and an output 

interface such as a screen, etc.) or any other electronic 
central unit, serving to send ultrasound compression 
waves into the medium 2 from its outside surface 3, which 
waves interact with diffusing particles 5 contained in 

10 the medium 1, which particles are reflective for 

ultrasound compression waves. The particles 5 may be 
constituted by any non-uniformity in the medium 1, and in 
particular, in a medical application, they may be 
constituted by particles of collagen present in human 

15 tissues (in echographic images, such particles form 
points known as "speckle") . 

To observe the propagation of the shear wave, an 
liltrasound probe 6 is used that is disposed against the 
outside surface 3 of the observed medium 1. This probe. 

20 delivers ultrasound compression wave pulses along an axis 
X, which pulses are of the type commonly used in 
echography, for example having a frequency lying in the 
range 0.5 MHz to 100 MHz, and preferably in the range 
0-5 MHz to 15 MHz, e.g. being about 4 MHz. 

25 The ultrasound probe 6 is constituted by an array of 

n ultrasound transducers Tl, T2 , Ti, Tn, where n_ is 

an integer not less than 1. 

By way of example, the probe 6 may be in the form of 
a linear strip capable of comprising, for example, n=12B 

30 transducers in alignment along an axis Y that is 

perpendicular to the axis X. However, the probe in 
question could equally be in the form of a two- 
dimensional array of transducers (plane or otherwise) . 
The transducers Tl, T2 , Tn are controlled 

35 independently of one another by the microcomputer 4, 

possibly via a central unit CPU which is contained for 
example in an electronics rack 7 connected via a flexible 
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cable to the probe 6. The transducers Tl-Tn can thus 
emit selectively: 

- either an ultrasound compression wave that is 
"plane" (i.e. a wave whose wave front is rectilinear in 

5 the X,Y plane), or any other type of unfocused wave 

illuminating the entire observation field in the medium 
2, for example a wave generated by causing random sound 
signals to be emitted by the various transducers Tl-Tn; 

- or else an ultrasound compression wave that is 
10 focused on one or more points of the medium 2 . 

To observe the propagation of the shear wave in the 
medium 2, several steps are performed in succession: 

a) an excitation step during which the microcomputer 
4 causes an elastic shear wave to be generated in the 

15 viscoelastic medium 2 by causing at least one ultrasound 
wave that is focused in the viscoelastic medium to be 
emitted by the probe 6; 

b) an observation step during which the propagation 
of the shear wave is observed simultaneously at a 

20 multitude of observation field points in the viscoelastic 
medium 2, this observation step comprising the following 
substeps : 

bl) the microcomputer 4 causes the probe 6 to 
emit into the viscoelastic medium a succession of 

25 unfocused ultrasound compression wave shocks at a rate of 
at least 500 shots per second (the focusing and the 
timing of the focus ultrasound wave emitted in step a) , 
and the timing of said unfocused ultrasound wave are 
adapted so that at least some of said unfocused 

30 ultrasound waves reach the observation field during the 
propagation of the shear wave through the observation 
field, for at least some of the unfocused ultrasound wave 
emissions ; 

b2) the microcomputer 4 causes the probe 6 to 
35 detect and record in real time sound signals received 

from the viscoelastic medium 2, said signals comprising 
echoes generated by the unfocused ultrasound compression 
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wave interacting with the reflecting particles 5 in the 
viscoelastic medium, these echoes corresponding (directly 
or indirectly) to successive images of the displacement 
of the viscoelastic medium; 
5 c) and at least one processing step during which: 

cl) the microcomputer 4 processes the 
successive sound signals received from the viscoelastic 
medium 2 during substep b2) in order to determine 
successive propagation images; and 

10 c2) the microcomputer 4 determines at least one 

movement parameter for the viscoelastic medium 2 at 
various points in the observation field. 

The focused ultrasound wave emitted during the 
excitation step a) may be a monochromatic wave of 

15 frequency f^ lying in the range 0.5 MHz to 15 MHz, for 

example being equal to about 4 MHz, which is emitted for 
a duration of k/ f seconds, where k_ is an integer lying in 
the range 50 to 5000 (e.g. being about 500) and f_ is 
expressed in Hz. Such a wave may possibly be emitted 

20 during a succession of emission periods separated by rest 
periods, the emission periods following one another at a 
rate lying in the range 10 to 1000 emissions per second. 

In a variant, the focused ultrasound wave emitted 
during excitation step a) is a linear combination (in 

25 particular a sum) of two monochromatic signals of 

respective frequencies fl and f2 such that 20 Hz < |fl - 
f 2 I < 1000 Hz, thus producing an amplitude modulated wave 
having a modulation frequency | fl - f 2 1 . 

In addition, the focused ultrasound wave emitted 

30 during excitation step a) may optionally be focused 

simultaneously or otherwise on a plurality of points so 
that the shear wave as generated presents a desired wave 
shape (for example it is thus possible to generate a 
shear wave that is plane, or on the contrary a shear wave 

35 that is focused) and illuminates desired zones in the . 
medium 2 . 
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During step bl) , which may last for example for less 
than one second, it is possible to emit unfocused 
ultrasound compression waves at a rate lying in the range 
500 to 10,000 shots per second, and preferably in. the 
5 range 1000 to 5000 shots per second (with this rate being 
limited by the go-and-return travel time for the 
compression wave through the medium 2, i.e. by the 
thickness of the medium 2 in the direction X: it is 
necessary for all of the echoes that are generated by the 

10 compression wave to have been received by the probe 6 
before a new compression wave is sent) . 

Each unfocused ultrasound compression wave 
propagates through the medium 2 at a propagation speed 
that is much higher than that of shear waves (e.g. about 

15 1500 meters per second (m/s) in the human body) , and 
interacts with the reflecting particles 5, thereby 
generating echoes or other analogous disturbances in the 
signal that are known in themselves under the name 
"speckle noise" in the field of echography. 

20 The speckle noise is picked up by the transducers 

Tl, Tn during substep b2) , after each shot of an 

unfocused ultrasound compression wave. The signal sij (t) 
as picked up in this way by each transducer Ti after shot 
No. j_ is initially sampled at high frequency (e.g. 3 0 MHz 

25 to 100 MHz) and is digitized in real time (e.g. on 

12 bits) by a sampler forming part of the rack 7 and 
connected to said transducer, the samplers being 
referenced respectively El, E2, En. 

The signal sij (t) as sampled and digitized in this 

30 way is then stored, likewise in real time, in a memory Mi 
belonging to the rack 7 and specific to the transducer 
Ti . 

By way of example, each memory Mi presents a 
capacity of about 128 megabytes (MB) , and contains all of 
35 the signals sij (t) received in succession for shots j = 1 
to p. 
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In deferred time, after all of the signals sij (t) 
corresponding to the same propagation of a shear wave 
have been stored, the central unit CPU causes these 
signals to be reprocessed by a summing circuit S 
5 belonging to the rack 7 (or else it performs this 
treatment itself, or indeed the treatment may be 
performed in the microcomputer 4) , using a conventional 
path- forming step corresponding to substep cl) . 

This generates signals Sj (x,y) each corresponding to 
10 the image of the observation field after shot No. j_ (when 
the unfocused ultrasound wave is a plane wave) . 

For example, it is possible to determine a signal 
Sj (t) by the following formula: 



15 where: 

- sij is the raw signal perceived by the transducer 
No. i after ultrasound compression wave shot No. j; 

- t(x,y) is the time taken by the ultrasound 
compression wave to reach the point of the observation 

20 field having coordinates (x,y) , with t = 0 at the 
beginning of shot No. j_; 

- di(x,y) is the distance between the point of the 
observation field having coordinates (x,y) and transducer 
No. i, or an approximation to said distance; 

25 - V is the mean propagation speed of ultrasound 

compression waves in the viscoelastic medium under 
observation; and 

- ai(x,y) is a weighting coefficient taking account 
of apodization relationships (in practice, in numerous 

30 cases, it is possible to assume that ai(x,y) = 1). 

The above formula applies muta^tis mutandis when the 
observation field is three-dimensional (with a two- 
dimensional array of transducers) , with space coordinates 
(x,y) being replaced by (x,y,2). 

3 5 After the optional path- forming step, the central 

unit CPU stores in a central memory M forming part of the 
rack 7, the image signals Sj (x,y) , or Sj (x) , or 



n 



Sj (t) 




ai(x,y) .sij [t (x,y) + di(x,y)/V] 



i=l 
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Sj(x,y,z), each corresponding to shot No. j_. These 
signals may also be stored in the microcomputer 4 if it 
is the microcomputer that performs the image processing 
itself. 

5 These images are then processed in deferred time in 

substep c2) by correlation and advantageously by cross - 
correlation either in pairs, or preferably with a 
reference image which may be: 

- either a displacement image determined previously 
10 as explained above and used as 'a reference image for 

subsequent displacement images (or for a limited number 
of subsequent displacement images), e.g. 30 displacement 
images) ; 

- or else determined during a preliminary initial 

15 observation step aO) , like the above-mentioned successive 
displacement images, by causing one or more unfocused 
ultrasound waves to be emitted by the probe 6 before 
excitation step a) which generates the shear wave (when a 
plurality of unfocused ultrasound compression waves are 

20 emitted in this way prior to the excitation stage, echoes 
generated by each unfocused compressed ultrasound wave 
are recorded interacting with the reflecting particles in 
the viscoelastic medium, these echoes corresponding to a 
plurality of successive preliminary images of the 

25 viscoelastic medium, and said initial image of the 
viscoelastic medium is determined by combining said 
successive preliminary images, and in particular by 
averaging the pixel values of said preliminary images) . 
The above-mentioned cross-correlation can be 

30 performed, for example, in a specialized digital signal 
processor (DSP) electronic circuit belonging to the rack 
7, or it may be programmed in the central unit CPU or in 
the microcomputer 4 . 

During this cross-correlation process, a cross- 

35 correlation function <S j (x, y) , S j +1 (x, y) > is maximized in 
order to determine the displacement to which each 
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particle 5 giving rise to an ultrasound echo has been 
subj ected . 

Examples of such cross-correlation calculations are 
given in the state of the art, in particular by O'Donnell 
5 et al. in "Internal displacement and strain imaging using 
speckle tracking", IEEE transactions on ultrasound, 
f erroelectrics, and frequency control. Vol. 41, No. 3, 
May 1994, pp. 314-325, and by Ophir et al . in 
"Elastography : a quantitative method for imaging the 
10 elasticity of biological tissues", Ultrasound Imag., 
Vol. 13, pp. 111-134, 1991. 

This produces a set of displacement vectors u(r',t) 
generated by the shear waves in each position of the 

medium 2 under the effect of the shear wave (these 
15 displacement vectors may optionally be reduced to a 
single component in the example described herein) . 

This set of displacement vectors is stored in the 
memory M or in the microcomputer 4 and can be displayed,, 
for example, in particular by means of the screen 4a of 

2 0 the computer, in the form of a slow motion picture in 

which the values of the displacements are illustrated by 
an optical parameter such as a gray level or a color 
level . 

The propagation differences of the shear wave 
25 between zones having different characteristics in the 

medium 2 can thus be seen clearly, for example the zones 
may comprise healthy tissue and cancerous tissue in a 
medical application . 

The motion picture of shear wave propagation can 

3 0 also be superposed on a conventional echographic image, 

which can also be generated by the apparatus 1 described 
above . 

Furthermore, it is also possible to calculate not 
the displacements of each of the points in the observed 
35 medium 2, but the defbrmations of the medium 2, i.e. 
vectors whose components are the derivatives of the 
displacement vectors respectively relative to the space 
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variables (X and Y coordinates in the example described) . 
These deformation vectors can be used like the 
displacement vectors for clearly viewing the propagation 
of the shear wave in the form of a motion picture, and 
5 they also present the advantage of eliminating 

displacements of the probe 6 relative to the medium 2 
under observation . 

From the displacement or deformation fields, the 
microcomputer 4 can advantageously then proceed with a 

10 map-making step d) during which, on the basis of the way 
in which the movement parameter (displacement or 
deformation) varies over time in the field of observation 
X, Y (or X, Y, Z with a two-dimensional array of 
transducers) , it calculates at least one propagation 

15 parameter of the shear wave, either at certain points in 
the observation field as selected by the user acting on 
the microcomputer 4, or else throughout the observation 
field. 

The propagation parameter of the shear wave that is 
20 calculated during the map-making step is selected, for 

example, from amongst: the propagation speed Cg of shear 
waves, the shear modulus /z, or Young's modulus E=3/i, the 
attenuation a of the shear waves, the shear elasticity 
jLtl, the shear viscosity //2, or the mechanical relaxation 
25 time Xg of the tissues. 

For example, it is possible at various points in the 
observation field to calculate: 

- the value of the propagation speed Cg of the shear 
wave, thus giving information about the hardness of the 

30 tissues; 

- the value of the mechanical relaxation time x g of 
the tissues, which is characteristic of the local 
viscosity of the medium. 

To do this, the following propagation equation (1) 
35 is used, with the displacement u(T,t) generated by the 
shear waves at each position r of the medium satisfying 
this equation: 
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dt^ dt 
where p is the density of the tissues, Xg is the 

mechanical relaxation time of the tissues, and c^ is the 

propagation speed of the shear wave which is directly 

5 related to Young's modulus E of the tissue by the 

following relationship : 

(2) C3 = J— 
V3p 

In the Fourier domain, the above wave equation (1) 
can be written as follows: 

10 (3) a>^pU(r,co) = cl (1 + jcox) . AU ( r , co) 

where U(r ,co) is the Fourier transform of the 
displacement field u(r,t) measured at each point, and 
AU(r ,co) is the Fourier transform of the spatial 
Laplacian of the field u(r',t). Given that COX3 << 1, it is 

15 possible to use a simplified expression: 

(4) =0} p 

(5) r = — tan(T( ) ) 

CO A\J{r, CO) 

where V|/ (x) is the phase of the complex variable x- The 

functions U(r ,0)) and AU(^,co) are known for each point of 

20 the echographic image, so it is possible to measure 

Young's modulus and the mechanical relaxation time of the 

tissue at each point in the space, thereby drawing up a 

map of those two parameters . 

Since equations (4) and (5) are true at each 

25 frequency, the calculation of Cg and x g can advantageously 

be averaged over the entire band of frequencies carried 

by the shear wave, thereby greatly improving the quality 

of the mapping that is performed. For this purpose, it 

is possible to use the folJ.owing formulae: 

30 (6) = i Im'p ia, 

f!l.^ il tam"^*!:"' ) ) 

cOy - cOq CO AU (r ,0)) 

where cOq and co^^ are* the minim and maximum frequencies 

carried by the shear wave . 

The method of calculation would be the same when 

35 using not displacements but deformations in the observed 

medium 2 . 
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Furthermore, it is advantageous to establish a 
succession of several maps of the desired propagation 
parameters, e.g. Cg and x^, by generating successive 
different shear waves, e.g. obtained by emitting 
5 ultrasound compression waves focused successively on a 
plurality of points or having different wave shapes. It 
is then possible to combine the various maps that are 
obtained in this way, e.g. by averaging them, so as to 
obtain a combination map that is richer and more 
10 accurate. 



